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Summary 

Optimal two-treatment, p period crossover designs for binary responses are determined. The 
optimal designs are obtained by minimizing the variance of the treatment contrast estimator over 
all possible allocations of n subjects to 2 P possible treatment sequences. An appropriate logistic 
regression model is postulated and the within subject covariances are modeled through a working 
correlation matrix. The marginal mean of the binary responses are fitted using generalized esti¬ 
mating equations. The efficiencies of some crossover designs for p = 2, 3,4 periods are calculated. 
The effect of misspecified working correlation matrix on design efficiency is also studied. 

Some key words: Binary response; Generalized estimating equations; Logistic regression; Efficiency. 


1. Introduction 

Crossover trials, wherein every experimental subject is exposed to a sequence of treatments over 
different periods of time, have been applied in a variety of areas; see for example Jones & Kenward 
(2014) and Senn (2003) for real life examples. The problem of determining optimal designs for 
crossover trials has been studied quite extensively in recent years and we refer to Bose & Dey 
(2009) for a review of results on optimal crossover designs. Most of the available results on optimal 
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crossover designs relate to situations where the response variable is continuous. However, there 
are situations in practice where the response in a crossover trial is binary in nature. For example, 
consider a trial reported by Senn (2003, page 127) wherein it was desired to study the effect of 
two drugs on 24 children aged 7 to 13 suffering from exercise-induced asthma. The two treatments 
were, a single dose of 12 pg formoterol solution aerosol, treatment A and a single dose of 200 pg 
of salbutamol solution aerosol, treatment B. Each child was given both the treatments either in 
the order, AB or BA. The response variable was binary in nature taking value 1 if the drug was 
effective and 0 otherwise. An equal number of children were allocated to each treatment sequence, 
AB or BA. Is this the best design to be used in this situation? 

In clinical or pharmaceutical research, the outcome of interest is often binary in nature. While 
methods for analysing binary data arising from crossover trials are available in Jones &; Kenward 
(2014) and Senn (2003), the question of designing such studies in an optimal manner does not seem 
to have been addressed. Waterhouse et al. (2006) considered crossover designs for binary response, 
where the treatments were taken to be continuous in nature and no period effects were considered 
in the model. 

In this article, optimal crossover designs are studied when the response variable is binary. We 
consider crossover trials with two treatments and p periods. The proposed designs minimize the 
variance of the estimator of treatment contrast of direct effects over all possible allocation of the 
n subjects to the treatment sequences. In the logistic regression model considered, both direct 
effect as well as the carryover effect of each treatment are considered, wherein we assume that the 
carryover effect of a treatment lasts only to the next succeeding period. While analysing data from 
binary crossover trials, often it is assumed that all observations are mutually uncorrelated; however 
this is not a very realistic assumption. We therefore assume that the p observations from each 
subject are mutually correlated while the observations from different subjects are uncorrelated. 
The correlation between observations within subjects are modeled using a “working correlation 
structure”. Since the main interest is in estimating the treatment effects, we treat the subject 
effects as a nuisance parameter and use the generalized estimating equations of Liang & Zeger 
(1986) to estimate the marginal means. Though estimating equations were used earlier by Jones 
&; Kenward (2014) to analyse repeated measures data, their models did not include the carryover 
effect. 
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The variance of the treatment effect estimator depends on the model parameters. To address 
the issue of parameter dependence, local optimal designs are found for given values of the model 
parameters. For p = 2,3,4, we study the effect of two working correlation structures, equi-correlated 
and autoregressive (AR) on the designs chosen. We also look at the effect of misspecification of the 
covariance on the design efficiency. 

In §2, we define the crossover logistic model for a binary response and discuss the estimation of 
the crossover model using generalized estimating equations. In §3, results on optimal two-treatment 
designs for 2, 3 and 4 periods are given. 

2. The model and estimation 

Consider a crossover trial involving t treatments, n subjects and p periods. Suppose the response 
obtained from the jth subject is Yj = (Y\j, ■ ■ ■, Y p j)', where a prime denotes transposition. Instead 
of specifying a joint distribution of the repeated measurements we use a working generalized linear 
model (GLM) to describe the marginal distribution of Yij as in Liang & Zeger (1986), 

f(Vij) = expliyijfaj - bifaj) + cfafajip], 

For a binary random variable Yjj, faj = log ^ —, b(faj) = log[l T exp {faj}}, c(yij ) = 0, and the 

t pij 

scale parameter ip is 1 (Robinson & Khuri, 2003). The mean of Y %] is pij and variance Pij(l — Pij). 
In a crossover setup, we model the marginal mean pij as 

logit(nij) = rjij = p + fa + t^j) + pd(i-i,j)’, i = 1 ■ ■ ■ ,p, j = 1, ■ ■ ■ ,n, (1) 

where p is the overall mean, fa represents the effect of the zth period, t s is the direct effect due to 
treatment s and p s is the carryover effect due to treatment s, , s = 1 ,t. Throughout, l u is a 
«xl vector of all ones, I u is the identity matrix of order u and 0 a b is an a x b null matrix. Also, 
we write /3 = (fa,..., fa)', r = (ri, - - -, T t )' and p = (p 1 ,..., p t )'. 

In matrix form the linear predictor corresponding to the jth subject, rjj = ■ ■ ■, T)pj)', can 

be written as 

Vj = ( 2 ) 
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where 9 = (p, /?, r, p)'. The design matrix is Xj = [l p Pj Tj Fj \, where Pj = I p \ T = (T {,..., T^) 7 , 
where Tj is a p x t matrix with its (*, s)th entry equal to 1 if subject j receives the direct effect of 
the treatment s in the zth period and zero otherwise; F = (F[ ,..., F^)' , where Fj is a p x t matrix 
with its (i, s)th entry equal to 1 if subject j receives the carryover effect of the treatment s in the 
ith period and zero otherwise. 

The estimating equations of Liang and Zeger (1986) and Zeger et al. (1988) are used to estimate 
the regression coefficients and to obtain their variances. It is assumed that measurements from 
the same subject in the p periods are correlated while observations from different subjects are 
uncorrelated. 

The dependencies between repeated observations from a subject are modeled using a ’’working 
correlation” matrix R(a ) where a is a vector of length s. If R(a) is the true correlation matrix of 
Yj, then 

Cov[Yj\ = A 1 / 2 R(a)A 1 / 2 , (3) 

I/O i /o 

where Aj = diag{p\j(l — pij), ..., p p j( 1 — p p j))- Let Vj = A-' R{a)A- . 

For a repeated-measures model, Zeger et al. (1988, equation (3.1)) define the generalized 
estimating equations (GEE) to be 

j = 1 


where p,j = (pij,.. ■, Ppj)' ■ The asymptotic variance for the GEE estimator 9 (see Zeger et al. 
1988, equation (3.2)) is 


Var{9) = 


vrl dpj 


3 = 1 


d9 i 89 


-i 


( 4 ) 


if CoviYj ) = Vj. However, if the true correlation structure varies from the “working correlation” 
structure, then Var(9 ) is given by the sandwich formula (Zeger et al. 1988, equation (3.2)) 


Var{9) = 


Y' d ^3 y-1 

^ d9 j d9 

3 = 1 


-1 


± B X V p Cm{ y l)V -^ 

3 = 1 


Y' ^T'j ! dpj 

d9 j 89 

3 =1 


-1 


( 5 ) 


For the crossover model (1), the ith element of ^ is = xVpij{ 1 — ), where xV is the 


4 










ith row of Xj for i = 1 ..... p. 

For finding optimal crossover designs for the logistic model we use the approximate theory as 
in Laska & Meisner (1985) and Kushner (1997, 1998). For a review of results on optimal crossover 
designs using the approximate theory, we refer to Bose & Dey (2009, Chapter 4). Fixing the 
number of subjects to n and periods to p, we determine the proportion of subjects assigned to a 
particular treatment sequence. Each treatment sequence is of length p and a typical sequence can 
be written as ui = (ii,..., t p )', U € {1,..., t}. Let 11 be the set of all such sequences. We denote 
by n u the number of subjects assigned to sequence ui. Then, n = > 0. A design ( 

in approximate theory is specified by the set {p u ,uj € 11} where p u = n u /n, is the proportion of 
subjects assigned to treatment sequence ui. 

Note that the matrices X) and Fj depend only on the treatement sequence uj to which the jth 
subject is assigned, so Tj = T u , Fj = F u - see Lemma 4.2.1 in Bose & Dey (2009). This implies, 
Xj = X u . Since np u subjects are assigned to treatment oj, the variance of 9 is 


Var c (6) = 


2^ P" d6 ^ dd 

.Ld(zQ 


-i -1 r 


2_^ n P“-gf V u Cov(Y u )V u — 




2^ P“ g 9 go 


n -1 


( 6 ) 


As before, if the true correlation of Yj is equal to R{a) then 


Var ( (6 ) 


E 


T r — i d F 

npuj^rV. 


i -i 


89 


89 


( 7 ) 


In crossover trials the main interest usually lies in estimating the direct treatment effect con¬ 
trasts. Thus, instead of working with the full variance-covariance matrix of 9 we concentrate on 
Var[f) where, 

Var c (f) = EVar ( (9)E\ (8) 


where E is a t x m matrix given by [Oti, 0 tp, E,0tt] and m is the total number of parameters in 9. 

An optimal design is one which minimizes the variance of f over the set of all possible allocations 
of the n subjects to the 2 P treatment sequences. 

3. Two TREATMENT CROSSOVER TRIALS: RESULTS AND DISCUSSION 


The optimality and efficiency of two-treatment crossover designs for continuous responses have 
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been studied by various authors, including Kershner & Federer (1981), Laska & Meisner (1985), 
Matthews (1987) and Carriere & Huang (2000). With two treatments of interest, the problem 
simplifies to minimizing the variance of the treatment contrast n — 72 to obtain optimal crossover 
designs. Reparametrizing model (1) as in Laska & Meisner (1985), using r = (t\ — T 2)/2 and 
P=(pi~ P2)/2 we get 

logit(nij) = n + /3i + + p<Z> d {i-i),j), (9) 

where 4? 4 = 1,3>b = — 1 and &d(o,j) = 0 . 

For illustration we go back to the example in §1, where there are two treatments, A and B 
applied in two periods to each child. The design used involved the treatment sequences AB and 
BA, with equal allocation to each treatment sequence. Thus, the matrix X u depends on the 
treatment sequence u € f2 = {AB, BA}. If the treatment sequence, for example is u = {AB}, then 


= 


/ 


V 


110 10 
101-11 


In the following, we look at the performance of the design, {AB, BA}. Senn (2003) fitted a logistic 
model with no carryover effect to the data set and computed confidence intervals for the various 
components of 8. Using these intervals we investigate if the above two-period design is the best 
choice in the given situation. We also look at general situations for determining optimal designs 
when p = 2,3 or 4 for the two treatment case. 

3.1 Designs compared. The designs that we consider are the same as those discussed by Laska & 
Meisner (1985) and Carriere & Huang (2000) and are listed below for p = 2, 3 and 4: 

(i) p = 2 : 

Design 1: AB and BA; Design 2: AB, AA, BA and BB, with equal number of subjects assigned 
to each sequence. 

(ii) p = 3: 

d \: ABB and BAA, with optimal allocation to each treatment sequence; 

^ 2 : ABB and BAA; 

g? 3 : ABB, AAA, BAA and BBB; 

c? 4 : ABB, AAB, BAA and BBA; 
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d§\ ABB, ABA, BAA and BAB-, 
d§: AAA and BBB. 

In designs d 2 — d&, each treatment sequence is allocated equally. 

(iii) p = 4: 

I: AABB, BBAA, ABBA, BAAB, with optimal allocation to each treatment sequence; 
II: AABB, BBAA, ABBA and BAAB; 

III: AABB, BBAA-, 

IV: ABBA, BAAB. 

In designs II-IV, each treatment sequence is allocated equally. 

For evaluating and comparing the above designs we define an efficiency measure as 



( 10 ) 


where is the optimal crossover design obtained. 

3.2. Working correlation structures. We consider the uncorrelated, compound symmetric or, equi- 
correlated and the AR(1) structures for the correlation matrix R(a). Under the equi-correlated 
covariance structure, Rj = (1 — a)I p + aJ p . 

Under the AR(1) assumption, 


Rj = cJ* * I, i ^ i! 


3.3 Results. We begin with the design for the trial reported in §1. A model without the carryover 
effect is fitted to the binary data. The parameter estimates and confidence intervals for the pa¬ 
rameters are: p € [1.1573,5.0893]; r € [-5.1738,-1.4029]; (3 £ [—2.8390,0.4932], as listed in Senn 


(2003). A parameter space for 6 is set up by taking Cartesian product of the confidence intervals. 
Under uncorrelated error structure, the efficiency of both designs 1 and 2 are the same. For the 
equi-correlated error structure, we take a = 0.2, 0.4 for finding the best design. The performance of 
the design {AB, BA} is studied through the distribution of the efficiency values in the parameter 
space and compared with the design {AB, BA, AA, BB}. During comparisons 10,000 values are 
randomly selected from the parameter space and the efficiencies computed. It is found that design 
1 is substantially superior to design 2 in terms of median and minimum efficiencies. Furthermore, 
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the efficiency of design 1 is less affected by an increase in the value of a. 

We next determine optimal crossover designs for p = 2,3,4 under the correlation structures 
stated earlier. For logistic models, the variance of t depends on the parameter values. To address 
the issue of parameter dependence we use various parameter spaces, £>i — Bq. For an explanation 
of the different parameter spaces considered, see Table CD For instance, suppose p = 2 and the 
parameter space is B\. Then, p £ [—0.5,0.5], /% £ [—1,1], 02 £ [—1,1], r £ [—1.5,1.5], p £ [—1,1] 
and a = 0.2. For each value of p, the designs listed in subsection 3.1 are compared by choosing 
10, 000 values from each of the six parameter spaces B\ — Bq. The performance of these designs are 
examined by noting their minimum and median efficiencies. The length of the parameter intervals 
considered in B\ — £>3 are more when compared to those from £>4 — Bq. This enables us to study 
the effect of increasing parameter uncertainty on the chosen designs. 

3.3.1. Estimation of direct treatment effects in the presence of carryover effects. When the carryover 
effect is included in the model and p = 2 , we only consider design 2 for estimation of direct treatment 
effects, because, as noted by Bose & Dey (2015), design 1 does not permit the estimation of contrasts 
among direct effects. First, let us consider the uncorrelated error structure. All designs for p = 2,3 
and 4 are compared over parameter spaces B\ and £> 4 . It is seen that design 2 has high median 
efficiency (99%) for p = 2 (see Table EJ). 

For p = 3 the distribution of efficiencies of the designs are presented in Table ED From this 
table we note that design d/± performs the best followed closely by designs d\ and Design d§ is 
the worst. The performance of design d^ is not much affected by increasing lengths of parameter 
intervals. Under uncorrelated errors for p = 4, design I performs the best in terms of median and 
minimum efficiencies (see Table HJ). On the average design I allocates equal proportion of subjects 
to each treatment sequence. 

We next consider a compound symmetric correlation structure for the errors. All designs were 
compared over six parameter spaces, £>1 — Bq. For a model including carryover effects and p = 2, 
design 2 has 99% median efficiency (see Table [2]). Thus, as in the case of continuous responses 
(Laska & Meisner, 1985), design 2 performs well even under a logistic model. Design 2 is also not 
much affected by changes in the correlation parameter and in the lengths of the parameter intervals. 

The distribution of efficiencies of the designs for p = 3 are presented in Table [3] From the table 
we see that with respect to median efficiency, the design d\ with optimal allocation of subjects to 


each sequence is the best, very closely followed by designs c ?2 and d±. The median efficiencies of 
these three designs are also not affected by increasing a. The minimum efficiency values of design 
c ?4 are the highest. The effect of increasing a is more on the minimum efficiencies values of designs 
d\ and d 2 as compared to d$. Thus, on the basis of both median and minimum efficiencies, the 
performance of design ^4 is the best among the six designs compared. Design dg has the worst 
efficiency values and its performance is affected by increasing values of a. Note that under the 
continuous response case, di with equal allocation to sequences ABB and BAA is known to be the 
universally optimal design within the class of three-period designs for equi-correlation (Laska & 
Meisner, 1985). 

For p = 4, median efficiencies of all designs I-IV are at least 98%, and do not change with 
increasing a (see Table 2 ]). However, when designs are compared with respect to minimum efficien¬ 
cies, design I is the best closely followed by design II, unchanged with changes in the correlation. 
Both designs III and IV record lower minimum efficiency values, which decrease with increasing a. 
On the average design I allocates re/4 subjects to each sequence. For the normal linear model case, 
Laska &; Meisner (1985) showed that design II is optimal. 

Finally we consider the AR(1) correlation structure. For p = 3, designs d\,d 2 and d^ record at 
least 99% median efficiency values (see Table EJ). However, the minimum efficiencies of designs d\ 
and d 2 are lower than those of design d^ and are also affected by increasing correlation. Even for 
the AR(1) structure, we see that the best design is dj on the basis of both minimum and median 
efficiencies and it is also least affected by increasing a. Design d$, with equal allocation to sequences 
AAA and BBB is the worst among the designs compared. Laska & Meisner (1985) showed that 
for p = 3 in the continuous response case, design d 4 is optimal. However for all a more individuals 
were allocated to sequence AAB and its dual. 

Under the AR(1) structure, for p = 4, design I is the best design with highest median and 
minimum efficiencies. Designs III and IV do not perform well, their minimum efficiencies are lower 
and decrease with increasing a (see Table [4]). The average allocations of design I is studied, the 
results show that for low values of a , the sequence ABBA and its dual each gets 35% allocation 
while, sequence AABB and its dual get 15% allocation. As a increases to 0.4, 85% individuals 
are allocated to ABBA and its dual, while for a = 0.6, it is 92%. Thus, as a increases almost no 
subjects are allocated to sequence AABB and its dual. Our result matches with those of Laska & 
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Meisner (1985) and Matthews (1987) when the response is continuous. 

3.3.2. Estimation of carryover effects. Although direct treatment effects are generally of primary 
interest in crossover trials, we present some results on optimal designs for estimating carryover 
effects as well for the uncorrelated, equicorrelated and AR(1) structures. For p = 3,4, the parameter 
spaces considered are B\ — B%. 

In the uncorrelated error case for p = 3, we note that d 4 performs the best with 99.51% median 
and 96.23% minimum efficiencies. The median efficiencies of d± and d% are also at least 99%; 
however their minimum efficiencies (90%) are much lower than those under d 4 (96%). 

For the equicorrelated and AR(1) covariance structures, d\ and c ?2 have 99% median efficiencies 
and at least 90% minimum efficiency values. Also the performance of designs d\ and c ?2 do not 
change with increasing correlation. However, unlike in the estimation of direct treatment effect the 
performance of d 4 is not good and is affected by changes in a. Median efficiency of <£4 dcreases 
to 94% and minimum efficiency to 90% as a increases to 0.6. Design d& is the worst among the 
designs compared as before. On the average d\ allocates equal (n/2) subjects to sequence ABB 
and its dual. For the normal case under AR(1) covariance structure, Matthews (1987) showed that 
d 2 has 100% efficiency for a = 0.2,0.4, 0.6. 

When p = 4, under all three covariance structures, design I performs the best with respect to 
minimum and median efficiencies. Under uncorrelated and compound symmetric covariances, the 
average allocation of design I is n/4 to each sequence. However, under AR(1), as a increases more 
subjects get allocated to sequence AABB and its dual than sequence ABBA and its dual. In the 
normal response case under AR(1) covariance structure, Matthews (1987) showed that design with 
sequences AABB, ABBA and their duals has at least 85% efficiency for positive a. 

3.4. Effect of covariance misspecification 

So far we have assumed that the true correlation structure of the responses is equal to the 
working correlation structure. However, this may not be true in most cases. To see the effect 
of varying structures on the efficiency of the designs, we carry out a simple study We set p = 3 
and use parameter spaces £>5 and Bq. The working correlation structure is taken to be compound 
symmetric and the true correlation AR(1). The designs studied are d\,d 2 and d 4 . The distribution 
of efficiencies of the three designs when working and true structures are assumed to be equal are 
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already given in Table 3. 

For finding the efficiency of the designs when there is misspecification we use (6) in the variance 
formula and report the results in Table [5j We see that the performance of all designs are affected; 
however designs d\ and c ?2 suffer much more than d 4 . In case of a = 0.4, the median efficiencies 
of designs d\ and ^2 reduce to 97%, and to 94% when a is 0.6. The efficiency of design d 4 is at 
least 97% for a = 0.4, 0.6. The minimum efficiencies of designs d\ and d 2 are affected much more 
than design d 4 and also more when we consider a higher value of a. Thus, design d 4 appears to 
be the best design to use when we have carryover effects in the model as it also guards against 
misspecification in covariance. 

3.5. Concluding Remarks. Crossover designs for binary responses are compared for p = 2,3,4. 
Since these designs depend on the parameter values, intervals of the parameters are considered and 
local optimal designs are found in each case. The main results on the estimation of direct effects 
are summarized below. 

For p = 3 design d 4 is seen to be the best design for uncorrelated, equicorrelated and AR(1) 
covariance structures. For p = 4 design I is seen to be the most efficient design for uncorrelated, 
equicorrelated and AR(1) covariance structures. In the equicorrelated case, on average design I 
allocates equal number of subjects to the treatment sequences. However, in the AR(1) case, the 
average optimal allocation depends on the correlation parameter a. The results found in the logistic 
regression case for p = 2, 3,4 are very similar to available results in the continuous case. 
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Table 1: Parameter spaces 


Parameters 


P 

Pi (for each i = 1, 2, 3, 4) 

T 

P 

a 


Bi 

[-0.5, 0.5] 

[-1.0, 1.0] 

[-1.5, 1.5] 

[-1.0, 1.0] 

0.2 


b 2 

[-0.5, 0.5] 

[-1.0, 1.0] 

[-1.5, 1.5] 

[-1.0, 1.0] 

0.4 


b 3 

[-0.5, 0.5] 

[-1.0, 1.0] 

[-1.5, 1.5] 

[-1.0, 1.0] 

0.6 


Bi 

[-0.5, 0.5] 

[-1.0, 1.0] 

[-0.2, 1.5] 

[-0.2, 1.0] 

0.2 


b 5 

[-0.5, 0.5] 

[-1.0, 1.0] 

[-0.2, 1.5] 

[-0.2, 1.0] 

0.4 


Bq 

[-0.5, 0.5] 

[-1.0, 1.0] 

[-0.2, 1.5] 

[-0.2, 1.0] 

0.6 


Table 2: Efficiencies (minimum, median) of design ( D 2 ) compared to the optimal design D* for 2- 
periods 2-treatments with carry-over effect model. Correlation Structures (Corr. st.) are compound 
symmetric (CS) and uncorrelated (Ind.). 


Par. Space Corr. St. D 2 


By 

CS 

(0.9639, 0.9977) 

Ind. 

(0.9599, 0.9974) 

b 2 

CS 

(0.9539, 0.9975) 

b 3 

CS 

(0.9595, 0.9972) 

Bi 

CS 

(0.9738, 0.9983) 

Ind. 

(0.9664, 0.9977) 

b 5 

CS 

(0.9749, 0.9981) 

Bq 

CS 

(0.9734, 0.9980) 
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Table 3: Estimation of direct treatment effect. Efficiencies (minimum, median) of designs d\ — d$ as compared to the optimal design D 
for 3-periods 2-treatnrents model. Correlation Structures (Corr. st.) are compound symmetric (CS), AR(1) and uncorrelated (Ind.). 


Par. Space Corr. St. (1\ 


d 2 


d 3 


G?4 


d§ 


d§ 



CS 

(0.9364, 0.9972) 

(0.9355, 0.9964) 

(0.9135, 0.9487) 

(0.9744, 0.9972) 

(0.9382, 0.9668) 

(0.7514, 0.8139) 

B\ 

AR(1) 

(0.9381, 0.9945) 

(0.9374, 0.9938) 

(0.9244, 0.9636) 

(0.9704, 0.9960) 

(0.9409, 0.9714) 

(0.7629, 0.8318) 


Ind. 

(0.9434, 0.9960) 

(0.9423, 0.9952) 

(0.9304, 0.9686) 

(0.9742, 0.9979) 

(0.9390, 0.9682) 

(0.7720, 0.8378) 

b 2 

CS 

(0.9293, 0.9974) 

(0.9291, 0.9966) 

(0.8964, 0.9308) 

(0.9707, 0.9959) 

(0.9337, 0.9658) 

(0.7255, 0.7796) 

AR(1) 

(0.9309, 0.9917) 

(0.9309, 0.9910) 

(0.9054, 0.9503) 

(0.9657, 0.9939) 

(0.9350, 0.9708) 

(0.7438, 0.8141) 

b 3 

CS 

(0.8955, 0.9971) 

(0.8955, 0.9963) 

(0.8830, 0.9156) 

(0.9659, 0.9943) 

(0.9247, 0.9650) 

(0.6803, 0.7315) 

AR(1) 

(0.9027, 0.9860) 

(0.9027, 0.9853) 

(0.8813, 0.9290) 

(0.9579, 0.9912) 

(0.9195, 0.9662) 

(0.6978, 0.7749) 

B a 

CS 

(0.9501, 0.9982) 

(0.9497, 0.9976) 

(0.9294, 0.9511) 

(0.9813, 0.9975) 

(0.9396, 0.9645) 

(0.7565, 0.8138) 

AR(1) 

(0.9467, 0.9961) 

(0.9461, 0.9954) 

(0.9417, 0.9666) 

(0.9821, 0.9970) 

(0.9415, 0.9690) 

(0.7680, 0.8328) 


Ind. 

(0.9501, 0.9973) 

(0.9489, 0.9966) 

(0.9446, 0.9725) 

(0.9853, 0.9986) 

(0.9403, 0.9651) 

(0.7699, 0.8375) 

b 5 

CS 

(0.9408, 0.9984) 

(0.9407, 0.9977) 

(0.9132, 0.9326) 

(0.9776, 0.9958) 

(0.9392, 0.9636) 

(0.7271, 0.7795) 

AR(1) 

(0.9310, 0.9937) 

(0.9308, 0.9931) 

(0.9233, 0.9533) 

(0.9773, 0.9953) 

(0.9369, 0.9689) 

(0.7478, 0.8160) 

Bq 

CS 

(0.9201, 0.9983) 

(0.9200, 0.9975) 

(0.8910, 0.9166) 

(0.9735, 0.9941) 

(0.9290, 0.9631) 

(0.6838, 0.7314) 

AR(1) 

(0.9036, 0.9881) 

(0.9034, 0.9875) 

(0.8894, 0.9317) 

(0.9672, 0.9929) 

(0.9158, 0.9647) 

(0.7016, 0.7777) 





Table 4: Estimation of direct treatment effect. Efficiencies (minimum, median) of designs I — IV 
compared to the optimal design D* for 4-periods 2-treatments model. Correlation Structures (Corr. 
st.) are compound symmetric (CS), AR(1) and uncorrelated (Ind.). 


Par. Space 

Corr. St. 

D t 

D n 

D ni 

Div 

Bi 

CS 

AR(1) 

Ind. 

(0.9803, 0.9996) 
(0.9772, 0.9994) 
(0.9767, 0.9985) 

(0.9787, 0.9982) 
(0.9746, 0.9962) 
(0.9759, 0.9970) 

(0.9151, 0.9848) 
(0.9046, 0.9728) 
(0.9180, 0.9840) 

(0.9228, 0.9845) 
(0.9315, 0.9927) 
(0.9245, 0.9838) 

b 2 

CS 

AR(1) 

(0.9805, 0.9996) 
(0.9754, 0.9999) 

(0.9792, 0.9983) 
(0.9678, 0.9906) 

(0.9000, 0.9845) 
(0.8816, 0.9571) 

(0.9164, 0.9842) 
(0.9204, 0.9972) 

b 3 

CS 

AR(1) 

(0.9733, 0.9996) 
(0.9720, 1.0000) 

(0.9707, 0.9982) 
(0.9592, 0.9829) 

(0.8894, 0.9837) 
(0.8591, 0.9404) 

(0.8860, 0.9837) 
(0.8975, 0.9986) 

B 4 

CS 

AR(1) 

Ind. 

(0.9878, 0.9997) 
(0.9870, 0.9996) 
(0.9875, 0.9989) 

(0.9847, 0.9986) 
(0.9837, 0.9973) 
(0.9855, 0.9977) 

(0.9504, 0.9876) 
(0.9369, 0.9761) 
(0.9520, 0.9872) 

(0.9186, 0.9811) 
(0.9275, 0.9911) 
(0.9228, 0.9812) 

b 5 

CS 

AR(1) 

(0.9917, 0.9997) 
(0.9807, 1.0000) 

(0.9840, 0.9987) 
(0.9770, 0.9926) 

(0.9368, 0.9872) 
(0.9163, 0.9607) 

(0.8999, 0.9810) 
(0.9151, 0.9969) 

Bq 

CS 

AR(1) 

(0.9899, 0.9997) 
(0.9720, 1.0000) 

(0.9834, 0.9988) 
(0.9644, 0.9852) 

(0.9306, 0.9868) 
(0.8835, 0.9448) 

(0.8859, 0.9809) 
(0.9030, 0.9988) 


Table 5: Efficiencies (minimum, median) of designs d\,d 2 and d 4 under misspecification of correla¬ 
tion structure. 

Par. Space 

di 

d 2 

di 

b 5 

(0.9124, 0.9716) 

(0.9120, 0.9710) 

(0.9651, 0.9821) 

B 6 

(0.8713, 0.9456) 

(0.8709, 0.9448) 

(0.9456, 0.9670) 
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